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We present polynomial force reconstruction from experimental intermodulation atomic force mi- 
croscopy (ImAFM) data. We study the tip-surface force during a slow surface approach and compare 
the results with amplitude-dependence force spectroscopy. Based on polynomial force reconstruc- 
tion we generate high-resolution surface property maps of polymer blend samples. The polynomial 
method is described as a special example of a more general approximataive force reconstruction, 
where the aim is to determine model parameters which best approximate the measured force spec- 
trum. This approximative approach is not limited to spectral data and we demonstrate how is can 
adapted to a force quadrature picture. 
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INTRODUCTION 

The combination of high-resolution imaging and high-accuracy force measurements [51llO| is a strong driving 
force for the development of atomic force microscopy (AFM) . The advent of multifrequency AFM resulted in a variety 
of new measurement techniques enabling enhanced contrast and spatial mapping of surface properties on a wide range 
of samples dlj. However, multifrequency AFM creates more data than conventional AFM, which both complicates 
the interpretation of measurement results and offers the possibility of much more detailed surface analysis. One of 
the goals when interpreting AFM data is the reconstruction of the force between a surface and the sharp tip at 
the end of the oscillating cantilever, while scanning. This reconstruction is readily possible by means of the Fourier 
transform if the tip's motion in response to this force and the cantilever's linear response function are known over a 
wide frequency band [T2Ul5j . However, the tip motion is often only measurable in a narrow frequency band around a 
cantilever resonance, since the cantilever transfer function sharply attenuates other frequency components of the tip 
motion, placing them below the detection noise floor. With this measured partial motion spectrum, the original force 
cannot be recovered with a simple Fourier transform and additional assumptions about the functional representation 
of the tip-surface force are required. 

These assumptions, which can be expressed with a finite set of parameters, result in a correlation of the measurable 
and the non-measurable frequency components of the motion. The parameters are chosen such that the spectrum of 
the reconstructed force best approximates the measured partial force spectrum. This approximative reconstruction 
requires the use of numerical optimization techniques if the force model is nonlinear in the parameters [TB]. Analytic 
solutions can be obtained if the model is linear in the parameters [TTHH]- Such a linear model of particular interest is 
the polynomial, as it constitutes a general expansion of the tip-surface force. 

Polynomial force reconstruction methods have been proposed theoretically and tested on simulated data for in- 
termodulation AFM fImAFM) jl71 [TO] . Here, we demonstrate, for the first time polynomial force reconstruction on 
experimental ImAFM data and compare it with reconstruction based on amplitude-dependence force spectroscopy 
(ADFS)[2D]. Moreover, fitting a force model to the polynomial reconstruction allows for the extraction of properties 
like surface adhesion, sample stiffness or interaction geometry. We demonstrate this extraction of surface properties 
with high resolution stiffness maps on a blend of polystyrene (PS) and poly (methyl methacrylate) (PMMA). 

Polynomial reconstruction, and most other multifrequency methods, work directly on the measured spectral data 
of the tip motion. Since the tip motion can be very complicated, the interpretation of spectral data often becomes 
rather difficult and alternative data representation schemes might provide a better understanding of the tip-surface 
force. Recently, we have shown how a narrow-band ImAFM measurement yields the oscillation amplitude dependence 
of a force component Fi in-phase with the sinusoidal tip motion and a force component Fq quadrature, or 90 degrees 
phase-shifted, to the tip motionl2II. Here we show how polynomial force reconstruction can be perfomed within the 
context of this picture of two force quadratures. 
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RESULTS AND DISCUSSION 



Polynomial force reconstruction from spectral data 

In narrow band AFM the tip dynamics as a function of time z{t) is usually described by a harmonic oscillator |221l23j. 
subject to an external drive force and a time-dependent tip-surface force 

Z+'^Z + k,{z-h)= FdHvo W + Ftsit) (1) 

where the dot denotes differentiation with respect to time, cjq, Q and are the mode's resonance frequency, quality 
factor and spring constant respectively, and h is the static equilibrium position of the tip above the surface. One 
should note that the time-dependence of the tip-surface force Fts can be considered as an implicit time-dependence, 
since it is assumed that the tip-surface interaction depends on the instantaneous tip position z and velocity z which 
are functions of time. In Fourier space equation ([T]) becomes 

Z{w) = X(W) (Fdrivc(w) + As(c^)) (2) 

where the linear response function 

Xi-) = ^2_^ + ,^ (3) 

determines the tip response to a sinusoidal force applied at the frequency uj. The drive force can readily be determined 
from a measurement of the tip motion far away from the surface, Zhcdt), where the tip-surface force is zero, 

^drivc(t^) = X""^('^)^frco(w) (4) 

If the broad-band tip response close to the surface z{uj) is known, one can easily solve equation ^ for the spectrum 
of the tip surface force 

Fts{c^) = X^^i(^) (^(w) - Zfi.cc(w)) (5) 

With the inverse Fourier transform, the time-dependent force acting on the tip can be readily determined from 
equation (|5|. 

Since the result of an experiment is a vector z of time-discrete samples of the continuous signal z{t) during a time 
window of length T = 1/Aw, the Fourier transform can be expressed using a unitary matrix F , 

z = Zz (6) 

z = z (7) 

where a single underline denotes a vector and a double underline a matrix. In a real experiment only a partial 
motion spectrum z^ can be measured since the cantilever's linear transfer function x suppresses response far away 
from resonance . Mathematically, this can be expressed with a diagonal windowing matrix operator that sets all 
frequency components outside the resonant detection band to zero. The measured spectrum is then given by 

Im = E 1 (8) 

and equation ([s]) for the measured data becomes 

(^m - Ifrce.m) = ^ X^' il - Ifrec) = W F,, (9) 



Since W is not invertable we cannot determine the complete force spectrum F^g from equation ([9|) and thus the time- 
dependent force remains unknown. To reconstruct the complete force spectrum from the measured partial motion 
spectrum Zm we expand the tip-surface force into a finite series from a set of basis functions 



N 



Fts{z,z) = ^fn{z,z)gn (10) 
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Our assumption that the force can be expanded in this manner results in a correlation of the unknown frequency 
components of the force with the measurable components. A common choice for the functions /„ to model conservative 
forces are monomials [17] 



but also other basis functions like 




i(z-/i)("+i)/2 modd 



(11) 



(12) 



for the representation of position-dependent viscosites have been considered [TO]. For a measured tip motion the force 
vector can the be approximated as 



Fts^Hg 



(13) 



where the coupling matrix H is given by 



( 



V 



(14) 



in which the columns are formed by the vectors /p, /j^, • • ■ , evaluated at the measured discrete tip positions and 
velocities. Here, we assume that the measured, or windowed tip motion is a good approximation of the complete 

(15) 



tip motion z. In Fourier space equation ( 13 ) becomes 
where 



H = H_{Zm,Zn-,) = 



^ f_Q (— m ' ) = Ll (— m ' ) ' ' ' = / jy (— m ' — m ) 



V 



(16) 



The force matrix equation ([9| can then be written as 
We introduce 

H H 



and solve equation (171 for g such that 



9=M1k M^m - ^frccm) 



(17) 



(18) 



(19) 



where _ff denotes the pseudo-inverse of .If there exists an unique solution for the coefhcient vector g, the matrix 



H equals the exact inverse of _ff .If there is more than one solution for g, equation ( 19 ) computes the solution for 



which the the vector g has minimum length. If no solution to equation ( 17) exists the pseudo-inverse approximates 
the inverse in a least-square sense. 

I using the Fast Fourier Transform (FFT) algorithm, 
le expansion coefficients g of the the tip-surface force. 



The matrix can be rapidly computed from equation ( 16 



Therefore, equation ( 19 1 provides a efficient way to determineT 



However, special careshould be taken to avoid aliasing effects due to the finite sampling of the data. To increase the 



numerical stability of equation ( 19 ) it is advantageous to normalize z^^ and z^ such that the largest absolute value of 



any vector element is 1. This normalization can be interpreted as a pre-conditioning procedure for the matrix H_ 
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To further investigate which information about the tip-surface force can be extracted, we focus on the monomial 




expansion basis defined in equation (11) and the case of narrow band ImAFM where the windowing matrix is given 
by 

fci < i < kM (20) 
else 

with Sij being the Kronecker delta, fciAcj the lower frequency limit of the resonant detection band and ku^^^ the 
upper limit. In figure jlj we plot the absolute values of the components of the matrix for experimental data. One 
could imagine applying different windowing matrices when building iJ^ , for example one which is weighted by the 
signal-to-noise ratio at each frequency. The absolute value of the element of -ff^with the index (i, j) is a measure of 
how much the j"^ element gj of the expansion coefficient vector g , contributes to the force at the i"^ frequency iAw 
in the force spectrum F^^. It is apparent from figure [l] that only polynomial coefficients of odd order contribute to 
the force measured in the resonant detection band when two drives close to resonance are used in ImAFM. Thus, 



equation ( 19 ) only yields the odd coefficients in the polynomial force expansion and the resulting polynomial force is 
odd with respect to z — h. To determine the missing even coefficients we assume that the tip-surface force is zero for 
z > ^non-interacting- With this assumption we fit the even polynomial coefficients while keeping the odd coefficients 
constant. This reconstruction method has been extensively tested and its accuracy verified with simulated datafTO]. 
In the following we will show results for experimental data. 



Polynomial force reconstruction during slow surface approach 

To demonstarte the capabilities of the polynomial force reconstruction we perform an ImAFM approach measure- 
ment on a silicon oxide surface. In this measurement two drive frequencies close to resonance result in a beat-like tip 
motion, with rapid sinusoidal oscillations and a slowly varying amplitude. The AFM z-piezo moves solowly towards 
the surface, such that during one beat period the static tip height above the surface can be considered to be constant. 

Figure [2] shows one frame from a movie (supporting information file 1) visualizing the measurement. For four 
consecutive beats in the time domain (a) the corresponding amplitude spectrum around the first resonance is displayed 
in (b) where the components or partial spectrum used for force reconstruction are marked with red circles. The 
polynomial force reconstruction is plotted (yellow solid line) in figure ^ together with an ADFS reconstruction using 
the same data (red circles) . In figure [2ji the amplitudes of the tip motion at the lower (red) and the higher drive 
frequency (yellow) are shown as functions of the z-piezo extension, and the vertical blue line indicates the current 
z-piezo extension. 
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FIG. 2: One frame from a surface approach movie (link to movie) showing the very onset of repulsive forces. The beating 
waveform (a) has the intermodulation spectrum around resonance (b) where red circles are analyzed to reconstruct the tip- 
surface force (c) using both both polynomial (yellow) and ADFS (red) methods. The z-piezo extension Zpiozo = 11-8 nm is 
indicated by the blue vertical line in (d) which the displays the amplitudes at the two drive frequencies. The interaction is 
purely attractive in the blue shaded area of (d) , becoming repulsive in the red shaded area. 



Far away from the surface the tip does not experience any surface force and the motion spectrum exhibits response 
only at the driven frequencies (supporting information file 1). Consequently, the reconstructed force is zero. As the 
surface is approached the attractive force regime due to the van-der-Waals forces between the tip and the surface 
is reached. In this regime new frequency components appear in the motion spectrum, so-called intermodulation 
products. Note that in the time domain the distortion of the signal is barely visible. Both polynomial and ADFS 
reconstruction show a growing attractive interaction until a force minimum of -1.75 nN is reached at a piezo extension 
of 11.8 nm. At this point the tip experiences hard mechanical impacts on the sample surface near the beat maximum 
which are manifest in the sharp onset of repulsive force in the polynomial and ADFS reconstruction (figure [2|;). 

As the z-piezo further extends the tip indents deeper into the surface and experiences stronger repulsive forces as 
shown in figure |3] where one frame of the movie at a piezo extension of 18.7 nm is shown. In this repulsive regime the 
polynomial and the ADFS reconstruction agree very well. However, the force minimum is slightly sharper shape with 
the ADFS reconstruction, which is not constrained to be continuous in high-order derivatives, as for a polynomial. 
The repulsive force reaches its maximum of 31 nN at a z-piezo extension of 23.6. Moving closer to the surface the 
maximum force during one beat decreases until the oscillation vanishes. 

During the whole surface approach the polynomial and the ADFS force reconstruction agree very well, indicating 
that the polynomial reconstruction accurately reproduces the force. The shape of the reconstructed force is very stable 
during the entire approach for both polynomial and ADFS reconstruction which is a result of the high signal-to-noise 
ratio for the measured frequency components close to resonance. The stability of the reconstruction during approach 
gives us confidence in the method's ability to accurately reconstruct sharp features in the force curve, like the force 
minimum. 
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FIG. 3: One frame from a surface approach movie (link to movie) showing strongly repulsive forces. The beating waveform 
(a) has the intermodulation spectrum around resonance (b) where red circles are analyzed to reconstruct the tip-surface force 
(c) using both both polynomial (yellow) and ADFS (red) methods. The z-piezo extension Zpiozo = 18.7 nm is indicated by the 
blue vertical line in (d) which the displays the amplitudes at the two drive frequencies. The interaction is purely attractive in 
the blue shaded area of (d) , becoming repulsive in the red shaded area. 



Surface parameter mapping 

Material scientists are often interested in determining surface properties with high spatial resolution. ImAFM 
measurements can be preformed while scanning a sample surface [24l I25j. enabling polynomial force reconstruction 
in every point of an AFM image. A specific force model can then be fitted to the complete force curve or parts 
thereof, generating a map of the model parameters over the complete surface. To demonstrate this we scanned the 
surface of a PS/PMMA blend with ImAFM. To the repulsive part of the polynomial force reconstruction we fit a 
Derjaguin-MuUer-Toropov (DMT) force model [26] of the form 

-^rcp(^) — -^min ^(^ ^min) ^ 

where Zmin is the position of the force minimum i^min and e is the DMT stiffness factor which depends on the tip 
radius and the effective stiffness of the tip-surface system. On should be aware of the fact that macroscopic force 
models like the DMT model might not be applicable on the nanoscale[27j and that tip shape and surface topography 
lead to an interaction geometry which is different from the model geometry of a perfect sphere and a perfect flat [20 j. 
Moreover, the DMT model does not account for adhesive forces in the contact regime which we try to circumvent by 
using tips with small radii. While the DMT model provides sufficient insight into material properties, the extracted 
numerical values of the DMT parameters should not be expected to agree with values for the bulk material. 

In fLgure|4]a map of the DMT stiffness factor is shown. Even thought, the two polymers are very similar in stiffness 
at room temperatures[25], two domains of different stiffness are clearly visible in the stiffness factor map. The stiffer 
domains are PMMA-rich and 10 nm higher than the surrounding matrix which is PS-rich and is a factor of two 
softer than the PMMA-domains. Similar results on the same model polymer system have been obtained with higher 
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FIG. 4: A surface property map showing the DMT stiffness factor e [GPa w}^^ ] with a histogram of the stiffness factor. 



harmonics methods [13 [29l and ADFSlSOl. 



Polynomial reconstruction from force quadrature data 



Polynomial force reconstruction is an approximative reconstruction method applied to spectral data obtained from 
Fourier analysis of the tip motion. The general idea, to determine the parameters of a force model such that an 
experimental observable is best approximated, is not limited to spectral data. An alternative to the spectral represen- 
tation is a picture which which represents the data in terms of two force quadratures. The force quadrature Fi is the 
component of the (time-dependent) tip-surface force that is in phase with the sinusoidal tip motion, the quadrature 
Fq is the force component which is phase-shifted 90 degrees to the tip motion [30]. 



Fi{A,Cj,h) = ^ [ Ft,iAcos{idt) 
^ Jo 



+ h, ~ujA sin(cji)) cos{u!t)dt 



FQ{A,Ld,h) 



T 



Fts {A cos{ujt) + h, —LoA sin(tjt)) sin(a;t)dt 



(21) 
(22) 



where the force quadratures are functions of the oscillation amplitude A, the oscillation frequency Co and the static 
tip height h, all of which are constant during each oscillation cycle. However, here we consider only the amplitude 
dependence of Fj and Fq which can be rapidly measured with ImAFM using a single oscillation cycle analysis that 
is based on a separation of time scales [21]. 

The representation of the measurement result in terms of the force quadratures Fi and Fq has the advantage that 
they are directly connected to the tip-surface force and independent of the actual complicated multifrequency tip 
motion. With spectral data certain points on the tip-surface force curve will receive greater weight if the tip spends 
more time at these positions. On the Fi{A) and Fq{A) curves the weight at each amplitude can be controlled by 
design. Furthermore, distortions due to feedback artifacts can easily be removed from the Fi{A) and Fq{A) curves 
and both conservative and dissipative forces can be analyzed separately. 

To demonstrate approximative force reconstruction on force quadrature data, we consider again a conservative 



polynomial force representation as in equation (11). For such a force Fq{A) — and equation (21 1 becomes 



1 ^ 

Fi{A) = -Y,g^A^^ 

n=0 



„n+l 



{ujt)dt 



(23) 
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FIG. 5: Comparison of force reconstructions for a tip interacting with silicon oxide. Polynomial reconstruction from spectral 
data (yellow) and force quadrature data (blue) and ADFS reconstruction (red circles) are in excellent agreement. 



The integral is non-zero only for odd n and by using 



cos"(0) 



n/2 



— y 



2"-i ^ V k 

k=0 



where 



is the binomial coefficient, equation (23) becomes 



N 



n = 
n odd 



9n 



A" 



cos((n- 2fc)0) 



N 

E 

n odd 



which implies that the odd polynomial coefficients of the force expansion can be obtained by simply rescaling of the 
coefficients of a polynomial approximation of the Fi{A) curve, 



1 



2"+i 



1 



To obtain the even coefficients we apply the same algorithm as for reconstruction from spectral data. 

We obtain a polynomial approximation of Fi{A) by simple polynomial fit to equidistant discrete points on the 
Fi[A) curve. Alternative methods like expansion in orthogonal polynomials or different types of interpolation poly- 
nomials with different convergence properties can also be applied. The polynomial reconstruction based on force 
quadrature data can be implemented even more efficiently than the rec onstruction on spectral data since multiple 
Fourier transforms to construct the coupling matrix Si in equation ( 16 1 are not required. 

To validate the equivalence of polynomial force reconstruction on spectral and force quadrature data, we consider 
the ImAFM approach measurement on silicon oxide described above. From the tip motion at a z-piezo extension of 
25.6 nm we compute the Fi{A) curve and remove all data points for which the oscillation amplitude was decreasing. 
From the polynomial approximation of the Fj[A) curve we obtain the force polynomial as described above. The 
resulting reconstruction is shown in figure [s] (blue line) together with the reconstruction from spectral data (yellow 
line) and an ADFS reconstruction (red circles). Over the full range of oscillation both curves agree extremely well. 
Excellent agreement is also observed at all other z-piezo extensions. 
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CONCLUSIONS 

Multifrequency AFM opens the window to a wide variety of novel AFM measurement techniques enabhng a much 
improved understanding of the forces between the tip and the surface. We showed that polynomial force reconstruction 
is an intuitive and powerful method to approximate this interaction, we demonstrated the method's use for accurate 
and detailed force measurement and for high resolution surface parameter mapping with experimental data. 

As the field of multifrequency AFM continues to evolve, new alternative data representation schemes can help to 
simplify analysis and extract more specific property of the tip-surface interaction. The force quadrature picture is 
such a scheme which decouples information about the tip-surface interaction from the actual tip motion. We showed 
how the general idea behind approximative force reconstruction can be adapted to the force quadrature picture and 
we introduced an efficient way to extract the polynomial coefficients from the force quadratures. 

We hope that in the future polynomial force reconstruction will be a useful method for many scientists and that 
new data representation schemes will inspire innovative analysis methods. 

EXPERIMENTAL 

The silicon oxide sample was cleaned in an oxygen plasma before performing measurements in a Bruker Dimension 
3100 AFM system. The cantilever (Bruker MPP-11120) was calibrated using a non-invasive thermal method[3T] and 
had a resonance frequency of fo = 229.802 kHz, a quality factor of Q = 396.9 and a spring constant of kc = 16.0 nm~^. 
The slow surface approach velocity was 2 nms^^. 

PS {M„ = 280 kDa, Sigma-Aldrich) and PMMA (M^ = 120 kDa, Sigma- Aldrich) were spin-cast from toluene 
solution with a concentration of 0.53 %wt at ratio of 3:1 (PMMA:PS). The sample was scanned in a Bruker Multimode 
2 AFM system with a cantilever BS 300A1-G (Budget Sensors) having a resonance frequency fo = 343.379 kHz, quality 
factor Q = 556.9 and spring constant kc = 35.1 nm~^. The maximum free oscillation amplitude close to the surfacce 
was 30 nm and we scanned an image with 256 x 1024 pixel within 17 minutes. 

For all measurements we used a intermodulation lockin analyzer (IMP 2-32, Intermodulation Products AB) which 
synchronizes the signal generation and acquisition for measurement of the multifrequency response [32J. 

SUPPORTING INFORMATION 

Movie showing the tip motion and the reconstructed tip-surface during a ImAFM approach measurement 
• Supporting information file 1 

— File name: sl.mp4 

— File format: mp4 

— Title: Tip-surface force during ImAFM approach measurement 
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